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■ The finite temperature Casimir free energy is caiculated for a dielec- 

trie ball of radius a embedded in an infinite medium. The condition 
e/i = 1 is assumed for the inside/outside regions. Both the Green func- 
tion method and the mode summation method are considered, and found 
to be equivalent. For a dilute medium we find, assuming a simple "square" 
dispersion relation with an abrupt cutoff at imaginary frequency il) = u>o, 
the high temperature Casimir free energy to be negative and proportional 
|_ to xo = Uoa. Also, a physically more realistic dispersion relation involv- 

ing spatial dispersion is considered, and is shown to lead to comparable 
' results. 



1 Introduction 

The Casimir effect has in recent years attracted considerable interest. For 
reviews on this topic see, for example, Plunien et aZ.Q], Mostepanenko and 
TrunovQ, and Milton ||. The problem has been rather difficult to solve, once 
curved boundaries are present, because of divergences in the formalism and in- 
sufficient contact with experiments. In this paper we will consider a compact 
ball of radius a, permittivity ei and permeability n±, situated in an infinite 
medium whose corresponding material constants are €2 and 
We shall take 

eiMi = £ 2M2 = 1 (1) 

so that the photon velocity is everywhere equal to c. 

The relativistic condition (1) was, in the context of the Casimir effect, to our 
knowledge first introduced by Brevik and Kolbenstvedt j|, ||. The reason why 
this particular condition was imposed was that when calculating the Casimir 
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surface force density on the sphere, in this paper called /, one avoids thereby the 
subtraction of a contact term: the contact term turns out to be zero just as in the 
case of a singular spherical conducting sphere in a vacuum Later, this kind 
of medium has been considered from various viewpoints [7 - 19]. A common 
feature of most of these papers is the use of macroscopic electrodynamics of 
continuous media. 

The main purpose of the present paper is to consider the Casimir effect 
for a compact ball satisfying the condition (1), at finite temperatures. Finite 
temperature Casimir theory has been considered previously by several authors, 
e.g., Lifshitz @ Fierz [fl), Sauer J§1, Mehra lljLBrown and Maclay pif, 
Milton et al. |25fl, Revzen et al. p6\ , and Klich [[19| . Also the very recent 
extensive treatise of Feinberg et al. |27| ought to be mentioned, since it describes 
in detail the characteristic properties of the high temperature limit, and contains 
also an extensive reference list. 

We start in the next section by considering, as a first step, the T = 
theory, and demonstrate the equivalence between the Green function method 
and the mode summation method. The equivalence has of course to hold if the 
theory is to make sense physically, but the equivalence is not so easy to see by 
a mere inspection of the mathematics. A detailed calculation is desirable. In 
section 3 we generalize the Casimir energy expression obtained in Jl5| to the 
fourth order in the diluteness parameter £; cf. Eq. (15) and the definition in 
Eq. (7). In section 4 we consider the Casimir free energy F(T) at finite T, 
emphasizing the difference between F(T) and the Casimir energy E(T), and 
give an order of magnitude estimate of F(T) in the case of a dilute medium 
endowed with a "square" dispersion relation implying a sharp high frequence 
cutoff at imaginary frequency Cj = ujq. The presence of a dispersion-induced 
strong, negative, cutoff-dependent contribution to the Casimir energy (or free 
energy, for finite T), seems first to have been pointed out by Candelas p8[ , for 
the case of a perfectly conducting shell. 

At high temperatures we find F(T) to be negative, decreasing linearly with 
T. This agrees with results found previously in [^6| and The Casimir 

energy E(T) itself tends to zero at high temperatures. 

2 Equivalence between the Green function ap- 
proach and the mode sum approach when T = 
0. Nondispersive medium 

In this section we consider first the Casimir surface force density / at zero 
temperature, for a nondispersive medium. There are in principle two different 
methods for calculating / (more precisely, the part of / that is independent of 
the high frequency cutoff): 

The first method is to make use of the Green functions. There are two scalar 
Green functions, corresponding to two different electromagnetic modes. Again, 
this method can be carried out in one of two different ways: 

1. First, one can evaluate / as the difference between the radial Maxwell 
stress tensor components on the outside and inside of the surface r = a. 
We have to subtract off the volume terms in the scalar Green functions, 
so that only the boundary dependent terms are left. This regularization 
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ensures that the surface force goes to zero when a ~ * oo, which is as 
we expect on physical grounds. Also, this regularization permits us to 
calculate the Casimir energy E once / is known, by the formula dE/da = 
—4ira 2 f. (At T = 0, E is the same as the free energy F.) 

2. Alternatively, one can calculate the energy E by integrating the Maxwell 
energy density over the volume, retaining, as above, only the boundary 
dependent terms in the Green functions. The outcome of these two ap- 
proaches are of course in agreement. 

The Green function method was first employed by Milton, for the case of an 
ordinary dielectric ball [^9|, and was later applied to the case of an e/i = 1 ball 
in Ref. {§. 

The result of the calculation can be given as follows. After a complex fre- 
quency rotation u — > ilu we can express all physical quantities as integrals 
over frequencies uj along the imaginary frequency axis.[] The cutoff dependence 
is omitted (it disappears in the regularization procedure). We define Riccati- 
Bessel functions si and ei by (x = ka) 

= J^Iu(x), ei(x) = \f^-K v (x), (2) 

corresponding to the Wronskian W{si,ei} = — 1. Here v = Z+ |; l v and K v 
are modified Bessel functions. We introduce two dispersion functions A™and 
A™ corresponding to the TE and TM modes: 

A™(x) = ^ei^ 2 sj(x)ei(x) - y^eiMi si{x)e' l (x), (3) 
Af M (x) = y/e^^(x)ei{x) - V^s;(x)e;(x). (4) 



The Casimir energy as derived in |Q on the basis of the Green function can 
be written as 

Here y = Lua, x = \y\\ the time difference t = t — t' has been Euclidean rotated, 
t — ► if, and 5 = f/a is the time-splitting parameter. Moreover, we have defined 
the function A;(x) as \i(x) = [s;(x)e;(x)]' [In Ref Q we used the symbols Di(x) 
and Di(x); the relations between these quantities and those defined by (3) and 
(4) are Di(x) = -^/JI^Aj E (x) and Di(x) = -^I/TJA™^). ] 

Consider next the second method, which consists in summing over the modes 
by means of contour integration. This method was recently used in Ref. fl5j , 
and is in turn based on the technique developed by Lambiase and Nesterenko 

1 This frequency rotation is a bit more involved than it might seem at first sight, since the 
Hankcl functions imply that there are singularities in the lower half of the complex frequency 
plane. Physically, these singularities occur because we assume an infinite exterior region 
(thus, no large exterior sphere making all eigenfrequencies real). The legitimacy of the complex 
frequency rotation in the presencauf these complexities follows from a more careful derivation; 
cf. the recent discussion in Ref. fefj . 
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and Nesterenko and Pirozhenko (3^] . One gets (when we supply the same 
parameter 5 as above) 

E = -iJ2( 2l + i ) /~ fH^M^WArw]- (6) 

Introducing the parameter £ by 

C = flp = ^, (7) 

the expression (6) can alternatively be written as 

if(a + i)£4^ta[i- { 'VW]. («) 

The two methods thus lead to two expressions for the Casimir energy, given 
by Eqs. (||) and (||). The expressions, of course, have to be equal, but to see 
the equality is not quite trivial by mere inspection. The equality is tantamount 
to the equality 

[A? E (x) A™ (*)]' = - ( ^ 2)2 [A? (*)] ' . (9) 

Now it turns out that this equality is not so difficult to prove after all if we 
insert the expressions (3) and (4) on the left hand side and take into account 
the above mentioned Wronskian between si and e;. Actually, the expression (6) 
was also made use of in the Green- function calculation of Milton and Ng pa] . 

3 A remark on the dilute ball when T=0 

Before leaving the zero temperature theory let us make a brief remark on the 
nondispcrsivc dilute ball, meaning that 

Kl « i- (io) 



The most practical expression to make use of for the Casimir energy, is that 
of Eq.(|J). We expand the logarithm to fourth order in £: 

in(i - exf) = -c 2 a?(x) - l -e>>t + o(e), (n) 

and make the following steps: 

1. put 5 = 0; 

2. perform a partial integration in Eq.(8); 

3. interchange summation and integration signs. 
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Then, omitting the 0(£ 6 ) terms in (fn]), 



1=1 jo 1=1 j ° 



(12) 



which can be processed further using the uniform asymptotic expansion for s/ 
and e/ up to 0(l/i/ 4 ), in the same way as in Refs. [|15| [j. There is one single 
divergent term in Eq. (12), which can be regularized by means of the Riemann 
zeta function, ((s). The only formula needed in practice is 



5>- s = (2 s -l)C(s)-2 s 



1=1 



(13) 



Moreover, we use the summations 
as well as the integral 
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(B is the beta function), to get 
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(15) 



This expression generalizes that of equation (3.10) in up to the fourth 
order in £. Numerically, the expressions between square parentheses in (15) 
are respectively 1.05966 and 0.88880. The fourth order contribution to the 
Casimir energy is thus weak and negative, corresponding to an attractive force 
component. 

In connection with this calculation, we wish to emphasise two points. First, 
we may avoid the time-splitting parameter 5 from the beginning, using zeta 
function regularization (here the Riemann function) instead. Secondly, the con- 
vergence parameter s that was for mathematical reasons introduced in fig] , can 
simply be omitted. At least from a pragmatic point of view, the nondispersive 
medium becomes quite simple. 



4 Dispersive medium: Finite temperatures 

We now take into account dispersive properties of the two media, meaning that 
e = e(ui), [i — (j,(uj). We shall still assume that e(oj)fi(u>) — 1, so that the velocity 
of photons in either medium is equal to the velocity of light. After frequency 
rotation, along the imaginary frequency axis, e = e(icD), (J. = fiiiCS). 

Some care ought to be taken when considering the dispersive formalism, 
since frequency derivatives of the material constants are no longer ignorable 
straightaway. Let us first consider the T = total zero-point energy, which we 
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shall call Ej + jj, of the inner (I) and the outer (II) regions. This energy is simply 
given by Eq. (6), which we repeat here with a slightly generalized notation: 

£ = -^i> + 1 ) /^|e^ln[AT E ^,«)A™^,a)]. (16) 

The reason for this equality is that the two dispersion functions Aj E and Af M 7 
as given by the dispersive generalization of Eqs. (3) and (4), retain their physical 
meaning also in the presence of dispersion. Thus the TE and TM modes are still 
determined by Aj E = and Aj M — respectively. The expression (16) needs 
to be regularized: we subtract off the zero-point energy E un if orrn corresponding 
to the limit a — > oo; cf. |jl^|. As si(x) — \e x and e\{x) — e~ x for large x, we 
get from (3) and (4) 

Af E (iu>, a -> oo) = Aj M {iQ, a -» oo) = + sfim). (17) 



The Casimir energy E = Ei + n — E un i / orm becomes accordingly, at T = 

~4Aj E {iu>,a)A™(iu;,a) 



(18) 



Observing that 1 + A; — 2s',ei, 1 — A; = — 2s;eJ we find (when reverting to the 
notation Aj E ' TM (x) instead of Aj E ' TAI (i<l>,a)) that 

Af E (x)Af M (x) = \(^T 2 + V^T) 2 [1 - exf(x)}, (19) 
which finally means that, at T = 0, 

5E(«+i)£|«£'»[i-e>(.w(.)]. (20) 

It is rather remarkable that this simple formula continues to hold even in the 
presence of arbitrary frequency dispersion. Since dispersion implies that we 
introduce a physically based high frequency cutoff, we expect that there is no 
longer any need for a time splitting parameter. We have therefore put 5 = in 
the expression (20). We will henceforth restrict ourselves to the first ( i. e., the 
second order) term in the expansion (11). 

In this section we adopt the simplest imaginable dispersion relation (a " square" 
form): we let the inner medium (a ball) correspond to 

/ ■ - \ fp = const, u<uj 
MM = i, . (21) 

[1, w > uj , 

whereas the outer medium is a vacuum. It is to be noted here that the general 
thermodynamic law saying that the susceptibility has to decrease monotonically 
along the positive imaginary frequency axis p3[ , is satisfied by this curve for 
e(iu>) — 1/ fj,(iu>) if it is given a small negative slope for u> < u)q. We still assume 
dilute media, |£| <C 1, as above. 

Consider now the case of finite temperatures. Since the thermal Green func- 
tion is periodic in imaginary time, we can replace the Fourier transform in the 
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T = theory with Fourier series in imaginary time. The transition to finite 
temperature theory is accomplished by means of a discretization of the frequen- 
cies, 

uj — ► uj n = 2im/ (3, x — > x n — Cj n a, (22) 

where n is an integer and (3 — 1/T (we put ks = 1)- The rule for going from 
frequency integral to sum over Matsubara frequencies is 

dx^t^ , t = 2na/f3. (23) 

n=0 

Here t is a nondimcnsional temperature, and the prime on the summation sign 
in (23) means that the n = term is counted with half weight. 

Since x n can be written as nt, we can express the Casimir free energy at 
finite temperatures as 

F ( r ) = ^E^E' A ^) ( 24 ) 

1=1 n=0 

to 0(£ 2 ). The corresponding zero temperature expression, still to be called E, 
is given by the first term to the right in Eq. (12). 

The following delicate point ought to be noted. We know that the mean 
energy per mode of the electromagnetic field is tp(u,fi) — \lo coth(i/3w). The 
function (p(oj, (3) is to replace ^w, when we construct the general contour integral 
for the Casimir energy E(T) at finite temperatures. It is the function coth(|/3w) 
which is responsible for the occurrence of the Matsubara frequencies along the 
imaginary frequency axis. Now, our formalism implies morover a transition from 
the energy E(T) to the free energy F(T). This is accomplished by means of a 
partial integration, making use of the relations 



dtu 8(3 7 dp 

For this reason we have changed the energy symbol from E to F in Eq. (24). 
At T = 0, obviously F(T = 0) = E(T = 0) = E. We thank I. Klich (personal 
communication) for helpful comments regarding this point. 

Let us in this context also recall how the surface force density / is calculated. 
Generally, as mentioned in section 2, / can be found as the difference between 
the radial Maxwell stress tensor components on the outside and inside of the 
surface r = a. The stress components are in turn constructed from the two- 
point functions like (Ei(r)Ek(r')) when r and r' are close together, at a given 
temperature T. The appropriate energy function is found by integrating 47ra 2 / 
over a, from initial position a = oo, to the final position where the radius 
of the sphere terminates at the fixed value a. This process is to take place 
at constant temperature. The appropriate energy function is accordingly the 
free energy F(T). Thus we see that, whereas in force considerations at finite 
temperatures one is lead naturally to the free energy, while starting from energy 
considerations one has instead to make use of thermodynamical considerations 
to calculate F(T). 
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Mathematically, the case n = requires special attention. For this purpose 
we observe the approximate expression for s; and ej at low x: 



si{x) 



I> 



^ /x\ 



ez(x) 



r(i/) /acx (-"+?) 



(25) 



From this it follows that si(x)ei(x) — x/2v, Xi(x) — l/2v, implying that the 
contribution to (24) from n = is 



F(T,n = 0) 



1=1 



(26) 



We add the contribution from n > 1, restricting ourselves to the dominant 
term in the uniform asymptotic expansion of A/ p5[ : 



Thereby 



F(T) 



(l + x 2 /v 2 )- 3 / 2 



2v 



87ra ^ v 
i=i 



l + O 



i 



2v b 



^ UP- + nH 2 f 

n—1 K ' 



(27) 



(28) 



We have here truncated the sum over I at I = Iq in order to avoid divergences. 
It is easy to give an order-of-magnitude estimate of lo: Our dispersion relation 
(21) implies that photons having frequencies u> higher than loq do not "see" the 
sphere at all. A photon of limiting frequency luq just touching the surface of the 
sphere has an angular momentum equal to ljqo, — xq. This type of argument 
has repeatedly been used in Casimir calculations [^8[ [l(], |ll[. One might put 
Iq = cxq where c is a coefficient of order unity. However, in view if simplicity, 
and since we are able to give only an estimate of the order of magnitude anyhow, 
we put henceforth c = 1. As upper limit in the summation we thus take 



l = Xq. 



(29) 



Now taking into account the exact summation formula derived in Eq.(A3) 
in S ( with p an arbitrary constant): 



E 



' (n 2 + p 2 y 



we can write Eq. (28) as 

C2 x ° 



640^ 



i=i 



■Kl> 

coth — 



1 3tt 

2p^ + W P 

37T 2 1 7T 3 coth7rp 

16p 4 sinh 2 np 8p 3 sinh 2 wp ' 



TTv/t 



2 /irv\ 2 coth7w/£ 



t sinh Tiv/t 3 V t J sinh ■Kv/t 



(30) 



(31) 



8 



the first term on the right hand side of (30) drops out. 

This expression can be calculated, with xq as an input parameter, to show 
how F(T) varies with t. It is of interest first to examine the case of low tem- 
peratures. Let us assume that 

t«l, (32) 

implying that nv/t 1. Then the first term in (31) dominates, and can be 
replaced by unity. We obtain approximately 

F{t ^° ) = ~^ X °- (33) 

This negative expression establishes a constant low-temperature plateau for 
the Casimir free energy. 

It is of interest to compare this result with that obtained from the slightly 
more complicated though physically more realistic Lorentz (or Sellmeir) disper- 
sion relation. The latter, when taken only along the imaginary frequency axis, 
can be written as 

utiui) - 1 = - — ^" . 9 = - — ° , (34) 

where xo is the zero frequency magnetic susceptibility. (The damping term is 
omitted.) Here ui = ui n represents a "soft" cutoff. This case was considered in 
0, for the zero-temperature case. For a dilute medium, \\a\ -C 1, the result of 
the calculation in |Q can be written 

*c-°>-*~£(*-i)- < 35 > 

This formula assumes mathematically, strictly speaking, that xq 3> 1 in addi- 
tion to t <C 1. The first condition is however not very stringent; it turns out 
numerically that (35) holds to better than 1 percent accuracy even when xq is 
as low as about 4. We thus see that the dominant first term in Eq. (35) agrees 
qualitatively well with Eq. (33), the latter obtained on basis of the "square" 
dispersion relation (21). The latter expression is 3/4 of the Lorentz-based ex- 
pression E ~ /16a)xo. This is a satisfactory result, in view of the semi- 
quantitative agreement between the meaning of the parameters xq occuring in 
(16) and in (34). 

The second term in Eq. (35), S(nondisp) = (3/64)£ 2 , is the characteristic 
positive Casimir energy (corresponding to a repulsive force component) for a 
nondispersive dilute ball satisfying the condition e/i = 1, at zero temperature. 
Cf., for instance, Q. 

The other limiting case of interest is that of high temperatures. Let us 
assume that 

f»x 0l (36) 

implying that irv/t -C 1 for all values of I in Eq. (31). Making use of sinh z = z 
and cosh z = 1 for small z, we obtain as leading term 

*<«»■»> -HE; ~SsE;- W 

1=1 r 1=1 
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This expression is independent of a. The high - temperature Casimir free 
energy, still negative, thus decreases proportionally with t. This is qualitatively 
in agreement with [Efl and p7[ . 




Figure 1: Using the "square" dispersion relation (21): Nondimensional free 
energy F re d(t) as defined by equation (38), versus nondimensional temperature 
t = 2na/P when xq = 10. 

Figures 1-3 show how the nondimensional reduced Casimir free energy 
F red (t), defined by 

F re d{t) = 17K~A V ' ( 38 ) 

3^s /(64a) 

varies with t, in the cases when xq — {10, 100, 1000}. The low-temperature 
plateaus evidently lie at F re d(0) = — 1. These figures are given on a linear scale, 
thus exhibiting the linear behaviour predicted by Eq. (37) at high temperatures. 
It is of interest to check the accuracy of (37) in cases when the condition (36) 
is reasonably well, but not extremely well, satisfied. Let us choose xq — 10, t — 
100 as an example. From the expression 

8t x ° 1 

Fred{t > X ) = V - (39) 

6irxn — v 
' 1=1 

we obtain, taking into account that 1/^ = 2.361749, that F re( i(xo = 10, t = 
100) = -20.047149. The full series solution yields the number -20.047154. The 
approximate formula (39) is thus in this case extremely accurate. Choosing 
t = 10, we get from Eq. (39) that F red (x = 10, i = 10) = -2.005, whereas the 
series (31) yields -2.091. Even in this case, where the condition (36) obviously 
is broken, the error in the formula (39) is thus only about 4 per cent. Generally 
speaking, at high temperatures when the diluteness of the medium becomes 
more pronounced, the Casimir free energy becomes small, as we would expect. 
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Figure 2: Same as figure 1, but with xq = 100. 



In connection with the temperature dependence of the free energy the fol- 
lowing physical argument should be noticed. From the dispersion relations (21) 
or (34) we would expect that the influence from the dispersion effect becomes 
important when the temperature becomes so high that the most significant fre- 
quencies in the thermal radiation field are of the same order as ujq- And as most 
significant thermal frequency, it is natural to take the value ui — uj rn correspond- 
ing to the maximum of the blackbody distribution. From Wien's dispacement 
law we know that u) m — 2.8/(3. Thus, we expect that dispersion becomes signif- 
icant when lo ~ loq ~ 2.8//3 which, in view of the definition t = 2 / na/[3 1 means 
that t ~ 2xq. Now, in Figs. 1-3 the linear scale chosen for the temperature 
implies that the intervals covered on the ordinate axes are small, and the men- 
tioned effect becomes hidden (the "shoulder" in each diagram appears to lie at 
a lower value of t). The effect becomes however clear if we replace the linear 
temperature scale by a logarithmic scale. Figure 4 shows, as an example, such a 
diagram for the case xq = 100. Then the ordinate interval becomes considerably 
larger, and we see that the prediction t ~ 2xo for the position of the "shoulder 
turns out to be reasonable. (In the diagram we have drawn the line t = 2xq to 
emphasize the transitional region.) The same behaviour was previously found 
to hold in [g. 

We note in passing that the high temperature entropy S turns out to be 
constant: from Eq. (37) we get 




(40) 



i=i 
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Figure 3: Same as figure 1, but with xq = 1000. 



5 An improved dispersion relation 

Instead of the formal divergence encountered in the previous section and the 
necessity to truncate the sum over I in Eq. (28) at an upper limit I = Iq, it 
is of interest to consider a mathematically more involved but physically more 
correct dispersion relation, that leads to a finite result for the free energy even 
after summing over all I up to infinity. 

Let us consider the following dispersion relation, being essentially of the 
Lorentz form but augmented by a term in the denominator that describes spatial 
dispersion: 



Xo_ 

l + (x 2 + l 2 )/xl ' 



M(w,0-l = q | ,„ 2 , I2W 2 , (41) 



Xo being a constant. As in the previous section we assume that the ball of 
radius a is dilute ( |xo | "C 1), and that the exterior region is a vacuum. The 
spatial dispersion term I jx\ ensures convergence when summing over all I. 
Since £ = (1 — /i)/(l + fJ.) according to (]?]), we obtain to first order in xo, at 
finite temperatures, 



-Xo/2 
l + (n 2 t 2 + l 2 )/x 2 



£ — *• 6n — 1 , /„ 2Ht 2 , n\i-2- (42) 



Instead of the expression (28) we now get 

2i/ 5 



t 2 r 4 f °° 
y ' 87TO ^ 



E 



. (43) 



v{x 2 + I 2 ) 2 ^ (x 2 + l 2 + n 2 t 2 ) 2 (v 2 + n 2 t 2 f 
where we have introduced the symbol £o = Xo/2- We resolve the last term into 
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Figure 4: Same as figure 2, but with a logarithmic temperature scale. Vertical 
line represents t = 2x$. 



partial fractions, and introduce the symbol S(p; k): 

oo 



(44) 



Then, if p and q are arbitrary constants, we get 
1 3 



y 

^— 4 (n 2 + p 2 ) 2 (n 2 + a 2 ) 

71=1 v y 



(p2 _ q 2y 



S(p; 1) 



(p 2 — g 2 ) 3 



S(p; 2) 



o z 



r 5(g;3). 



(p 2 — g 2 ) 2 

From Appendix A in |gj we write down the expressions for 5*(p; fc); fc = 1, 2: 

1 7T 



(45) 



SO?! 1) = -TTo + ^-coth?rp, 
2p 



(46) 



„. 1 7T 7T 2 1 

b(p; 2) = -—7 + — - cotIi7rp + — 5 3 

2p A 4p 3 y 4p 2 s inh 2 7rp 



(47) 



(S(p;3) was given in Eq. (30)). We define, in analogy with Eq. (38), the 
nondimcnsional reduced free energy to be 



Fred(t) 



F(T) 



3#c„/(64o)' 



(48) 
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and obtain after some manipulations the following expression: 

1 



3tt ^\2v(xt + P) 
Zv 5 



t\xl-l-\f 

+ ^-V{F S( I ;3) } 



5(^±^;l)-^;l) 
5(^±?;2) + 2S(^;2) 



(49) 



This expression is calculated numerically for two of the cases above, viz. x 
{10, 100}. The results are shown in Fig. 5. 



-0.6 




Figure 5: Using the spatial-dispersion relation (41): Nondimensional free 
energy F re d{t) as defined by Eq. (48), versus nondimensional temperature 
t = 2Tta/[3 when x n — {10, 100}. The low-temperature plateau is in accor- 
dance with Eq. (54). 

It is seen that the curves have roughly the same form as in Figs. 1-2. 
There is a definite low - temperature plateau, and for high values of t the 
curves tend towards a linear form. In any case, whatever we adopt a square 
dispersion relation or a spatial - dispersion relation, the sensitivity with respect 
to increasing values of t is most pronounced when xo is smallest. The limiting 
values for t — + in Fig. 5 are in agreement with those of Figs. 1 and 2 to within 
about 70 per cent. This is the order of magnitude - agreement that we might 
expect, since the parameter £ made use of in Eq. (48) is similar to, but not 
identical with, the nondispersive quantity £ made use of in Eq. (38). 
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It is of interest to give analytic approximations for the case of low tempera- 
tures, Then the quantities p = y/ Xq + I 2 jt and q — v/t are large, and it 
becomes convenient to make use of the approximate formula [|[ [H| 

y 1 0Fr(fc-l) 1 fn\ k e-^ 

^{n 2 +p 2 ) k 2T{k)p 2k - 1 2p 2k \p) T{k) ' 1 ' 

from which it follows, in our notation, that 



Inserting Eqs. (51)-(53) into Eq. (49) we obtain a reasonable simple series 
expression for F rec i(t), at low temperatures. A further simplification can be 
obtained if we keep only the leading order terms in Eqs. (51)-(53). For p ^> 1 
this means that we keep only the first term in each equation. We then get in 
this limit 




ed(t -» 0) 



In this expression, the temperature t is no longer present. Numerical trials 
show that Eq. (54) is very accurate for t <C 1. Choosing xq = 10, t = 1, we 
obtain from Eq. (54) F re< i(xo = 10, t = 1) = —0.6526, whereas the full formula 
(49) yields the number -0.6530, i. e., an error of 0.06 per cent. Even with t = 2, 
where Eq. (49) yields the number -0.678, the error in Eq. (54) is increased to 
only about 3 per cent. 

6 Conclusion and final remarks 

Let us summarize as follows: 

(1) For a relativistic medium, i.e., a medium satisfying the condition e/i = 1, 
the T = result for the Casimir energy is most conveniently written in the form 
(8). This expression holds for a nondispersive medium, S — r/a being the time- 
splitting parameter. The permittivity is here arbitrary; the medium need not 
be dilute. The regularization is made by subtracting off the volume terms in the 
two scalar Green functions. This regularization is natural, among else things 
because it permits one to relate E to the surface force density / by the relation 
dE/da — —4ira 2 f. The equivalence between the Green function approach and 
the mode sum approach is demonstrated explicitly. 
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(2) As an additional by-result at T = 0, the expression (15) generalizes the 
second order result for E found in |l5[ up to the fourth order in the parameter 

(3) Including the dispersive effect, the simplest dispersion relation one can 
imagine is the "square" relation given in Eq. (21). The finite-temperature 
expression for the free energy on the basis of this dispersion relation is given 
in Eq. (31), with Xq = uj^a being a high frequency cutoff which is in practice 
dependent on the detailed structure of the medium. With reasonable accuracy 
xq can represent the maximum Iq of the angular momentum variable I. Typical 
results are shown in Figs. 1-4. In the limiting case t — > they are in accordance 
with Eq. (33), and for f > io they are in accordance with Eq. (39). In 
particular, the free energy is for high temperatures negative, and is a linear 
function of t. The entropy is correspondingly a (positive) constant. These 
results are in accordance with those of Refs. p^j and |2?[| . The transition region 
between low-temperature and high-temperature theory can be determined with 
reasonable accuracy from a simple physical argument. 

These results are similar to those obtained on the basis of the Lorentz dis- 
persion relation in || . 

(4) One can take spatial dispersion into account, as we have done in section 
5. The basic dispersion relation is then Eq. (41). A physically appealing feature 
of this kind of approach is that the sum over angular momenta does not have 
to be truncated. Typical results for this case are shown in Fig. 5. 

(5) The condition e/i = 1 simplifies the theory, as noted already in Ref. Q. 
For an ordinary dielectric ball (/i = 1) the calculation becomes more complex 
and difficult to interpret p5| . 

(6) A final remark is called for, as regards the numerical computation of 
the series in the case of spatial dispersion. The evaluation of the expression 
(49) turned out to be more difficult than one might expect beforehand. Thus, 
a simple use of Matlab turned out to be insufficient. The problems seem to 
be associated with numerical noise, caused by the "critical" terms for which 
I ~ x§. There are relatively large individual terms that almost, but not quite, 
compensate each other in the sum. We managed to do the calculation making 
use of double - precision MS-DOS Quick Basic. In practice, adequate precision 
was found with inclusion of less than i„ terms. 
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